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ABSTRACT 



i— H Aims. We present coordinated near-infrared AT-band interferometric and optical spectroscopic observations of the M 1.5 giant aCet (Menkar) 
\Q ■ obtained with the instruments VINCI and UVES at the Paranal Observatory. Spherically symmetric PHOENIX stellar model atmospheres are 
' constrained by comparison to our interferometric and spectroscopic data, and high-precision fundamental parameters of Menkar are obtained. 
Methods. Our high-precision VLTI/VINCI observations in the first and second lobes of the visibility function directly probe the model- 
Q | predicted strength of the limb darkening effect in the K-band and the stellar angular diameter. The high spectral resolution of UVES of 
I ' R = 80 000 - 1 10 000 allows us to confront in detail observed and model-predicted profiles of atomic lines and molecular bands. 

Results. We show that our derived PHOENIX model atmosphere for Menkar is consistent with both the measured strength of the limb-darkening 
in the near-infrared A"-band and the profiles of spectral bands around selected atomic lines and TiO bandheads from 370 nm to 1000 ran. At 
the detailed level of our high spectral resolution, however, noticeable discrepancies between observed and synthetic spectra exist. We obtain 
a high-precision Rosseland angular diameter of 0r oss = 12.20 mas ±0.04 mas. Together with the Hipparcos parallax of 14.82 mas ± 0.83 mas, 
^ ' it corresponds to a Rosseland radius of i?R OSS =89 ± 5 R G , and together with the bolometric flux based on available spectrophotometry, to an 
k> effective temperature of T e g =3795 K + 70 K. The luminosity based on these values is L= 1460 L Q ± 300 L Q . Relying on stellar evolutionary 
. tracks, these values correspond to a mass M = 2.3 M G ± 0.2 M Q and a surface gravity log g=0.9 ± 0. 1 (cgs). 
C\3 Conclusions. Our approach illustrates the power of combining interferometry and high-resolution spectroscopy to constrain and calibrate 
stellar model atmospheres. The simultaneous agreement of the model atmosphere with our interferometric and spectroscopic data increases 
confidence in the reliability of the modelling of this star, while discrepancies at the detailed level of the high resolution spectra can be used to 
further improve the underlying model. 

Key words. Techniques: interferometric - Stars: late-type — Stars: atmospheres - Stars: fundamental parameters - Stars: individual: a Cet 
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1. Introduction 

Stellar atmosphere models predict in general the spectrum 
emerging from every point of a stellar disc. These models 
are often based on assumptions of a stationary 1 -dimensional 
plane-parallel atmospheric structure (e.g., ATLAS 9 models, 
Kurucz 1993), or a stationary spherically symmetric structure 
(e.g., PHOENIX models, Hauschildt & Baron [T999l Hauschildt 
et al. rH%9l . 

Confrontation of model atmospheres with observations is 
often performed by spatial integration of the model radiation 
field over the whole stellar disc, spectral integration to broad- 
band colours and/or spectrograms of different spectral res- 



* Based on observations collected at the European Southern 
Observatory, Chile, under program number 71.D-0370(A,B). 



olution, and subsequent comparison to observed broadband 
photometry, spectro-photometry, or (high-resolution) spec- 
troscopy. For instance, Tripicchio et al. dl997l fT999l measured 
the D resonance lines of neutral sodium and the KI resonance 
line at 7699 A to test MARCS (Gustafsson et al. 1 19751 Plez et al. 
1992 1993 1 model atmospheres that include photospheric and 
chromospheric effects. Valenti et al. (fl998 1 performed a least- 
squares fit of synthetic spectra to echelle spectroscopy of an 
M3.5 dwarf. They successfully constrained the model atmo- 
sphere by using the 7087 A TiO bandhead and five atomic Fe I 
and Til lines between 8674 A and 8692 A. Decin et al. (2003 1 
compared infrared spectroscopy obtained with ISO-SWS with 
MARCS stellar model atmosphere predictions to estimate funda- 
mental stellar parameters. 
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Optical interferometry provides a further and more direct 
test of stellar atmosphere models by resolving the stellar disc 
and measuring the centre-to-limb intensity variation (CLV) 
across the stellar disc. Interferometric measurements of stel- 
lar CLVs include, for instance, those by Hanbury Brown et 
al. (fT974t . Haniff et al. (FT9951 . Quirrenbach et al. (fT996t . 
Burns et al. J1997i . Hajian et al. (11998b . Wittkowski et al. 
d200Tl 120041 l2006t . Perrinet al. (I2004al l2004bl . Woo druff e t 
al. Fedele et al. ( 120051 . and Aufdenberg et al. (12005V 

The combination of interferometry and high-resolution spec- 
troscopy to constrain stellar model atmospheres, as suggested 
and presented, for instance, by Aufdenberg et al. |2003l and 
van Belle et al. ( 2003 1, has so far been surprisingly rare, al- 
though this approach can provide a very strong test of theoreti- 
cal models. 

Here, we present coordinated near-infrared Tif-band inter- 
ferometry obtained with the VLTI instrument VINCI and high- 
resolution (R up to 1 10 000) optical spectroscopy obtained with 
the VLT echelle spectrograph UVES of the M 1.5 giant a Cet 
(Menkar, HD 18884, HR911, HIP 14135). We demonstrate 
the advantage of combined interferometry and high-resolution 
spectroscopy to test and constrain PHOENIX model atmospheres 
(Hauschildt & Baron fT999l Hauschildt et al. lT999l We derive 
a set of high-precision fundamental parameters of a Cet. 

a Cet is a red giant with spectral type M 1.5 III (Morgan 
& Keenan fl973> . a parallax of 7r=14.82mas ±0.83 mas, a 
V magnitude of 2.54 (both from Perryman & ESA I1997> . 
and a small variability of AV=0.09 (Kholopov et al. 1999 1. 
Decin et al. (I2003i derived an effective temperature of 
r eff =3740 + 140 K, a surface gravity of log g=0.95 + 0.25, 
a mass of M/M =2.69 + 1.61, a [Fe/H]= 0.00 + 0.30, a 
luminosity of L=1455 + 328 L , and a (limb-darkened) an- 
gular diameter of 12.52 + 0.79 mas, based on a com- 
parison of ISO SWS data to MARCS models. Alonso et 
al. J2000> derived a (limb-darkened) angular diameter of 
12.59 + 0.36 mas and an effective temperature of 3704 + 39 K 
by means of the infrared flux method (IRFM). Earlier in- 
terferometric diameter measurements include those by Dyck 
et al. J 19981 ©ud at K-band 11.6 + 0.4 mas; bolometric 
flux 1.05+0.16 10~ 12 WcirT 2 ), Quirrenbach et al. (tT993l UD 
at 712nm, TiO band, 1 1.95 +0.23 mas; 0ud at 754nm, 
continuum, 11.66 + 0.22 mas), Mozurkewich et al. ( 119911 
120031 UD at 800nm 12.269±0.237mas; UD at 550nm 
1 1 .473 +0.25 1 mas; UD at 45 1 nm 1 1 .325 +0.4 1 mas). 

2. VINCI observations and data reduction 

We obtained near-infrared A'-band interferometric measure- 
ments of a Cet in the first and second lobe of the visibility 
function using the ESO Very Large Telescope Interferometer 
(VLTI) equipped with the VINCI instrument and the two VLTI 
siderostats. For a recent technical description of the VLTI, see 
for instance Glindemann et al. (|2003 1; for the VINCI instru- 
ment, see Kervella et al. J2003ll2004i . 

Between 2 October 2003 and 5 December 2003, 38 visibil- 
ity measurements were secured in the first lobe of the visibility 
function using the baseline E0-G0 with an unprojected baseline 
length of 16 m. Further observations obtained during the same 



Table 1. Properties of the used interferometric calibration stars. 
The references are (D) calibration by Dyck et al. (19961, (B) 
Borde et al. (20021, and (K) Kervella et al. (120031 . 



Calibrator ©ud Ref. Sp. Type T t g 



17Mon 


2.59+0.26 


D 


K4III 


4090 


18Mon 


1.86+0.02 


B 


KOIIIa 


4656 


19Ari 


2.37+0.03 


B 


M0III 


3690 


28 Mon 


3.14+0.30 


D 


K4III 


4090 


29 0ri 


1.92+0.19 


D 


G8III 


4670 


30 Gem 


2.02+0.03 


B 


K1III 


4581 


31 Leo 


3.22+0.05 


B 


K4III 


4202 


310ri 


3.56+0.06 


B 


K5III 


4046 


6 Leo 


2.10+0.03 


B 


K2.5IIIb 


4318 


/3Cet 


5.18+0.06 


B 


K0III 


4656 


^Phe 


2.69+0.03 


B 


K5III 


4046 


HR 1663 


2.69+0.27 


D 


K5III 


3920 


7] Cet 


3.35+0.04 


B 


K1.5III 


4380 


j/Eri 


2.64+0.20 


D 


Kim 


4510 


y'Cae 


2.31+0.23 


D 


K3III 


4230 


a- Car 


6.45+0.60 


D 


F0II 


7200 


HR 1799 


2.04+0.20 


D 


K5III 


3920 


HR2311 


2.43+0.04 


B 


K5III 


4046 


HR4546 


2.53+0.04 


B 


K3III 


4256 


HR6862 


2.59+0.05 


B 


K4.5III 


4097 


HR7092 


2.81+0.03 


B 


M0III 


3690 


tCet 


3.27+0.04 


B 


K1.5III 


4508 


i Eri 


2.12+0.02 


B 


K0III 


4581 


l Hya 


3.41+0.05 


B 


K2.5III 


4318 


^CMa 


2.38+0.03 


B 


Kim 


4497 


ttCMa 


5.93+0.02 


K 


A1V 


9230 


0CMa 


4.04+0.40 


D 


K4III 


4099 


<THya 


3.33+0.30 


D 


G9II-III 


4559 



period (11 additional nights during October to December) and 
using the same 16 m baseline were not used because a consis- 
tent interferometric transfer function could not be established 
with sufficient precision, partly owing to an insufficient time 
coverage of calibration stars. On 11 & 12 January 2004, four 
visibility measurements were recorded in the second lobe of 
the visibility function employing the D0-H0 baseline with an 
unprojected length of 64 m. Since all our observations cover 
only a small range of projected baseline angles on the sky (see 
Table |3, they are not sensitive to possible asymmetries of the 
stellar disc. Data of a Cet and of several interferometric cal- 
ibration stars were obtained as series of 500 interferograms 
with scan lengths of 250/im and a fringe frequency of 295 Hz. 
Calibration star details, including the adopted angular diame- 
ters, are listed in Tabled 

Mean coherence factors were obtained for each series of 
interferograms using the VINCI data reduction software, ver- 
sion 3.0, as described by Kervella et al. (2004), employing the 
results based on the wavelet transforms. Calibrated squared vis- 
ibility values for a Cet were obtained from the mean coherence 
factors as described in Wittkowski et al. (2004i. A time kernel 
of 3 hours was used to convolve the measured interferomet- 
ric transfer function during each night. The computed errors of 
the squared visibility amplitude take into account the scatter of 
the single scan's coherence factors, the errors of the diameters 
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Table 2. Details of our VINCI observations of a Cet (date 
and time of observation, spatial frequency, azimuth angle of 
the projected baseline (E of N)), together with the measured 
squared visibility amplitudes and their errors. The last column 
denotes the number of successfully processed interferograms 
for each series. 



Date & 


Sp. freq az 


V 2 


cr V 2 


# 


Time (UT) 


[1/"] [deg] 









of the calibration stars, and the variation of the interferometric 
transfer function over the night. Observational details and the 
calibrated squared visibility values with their errors are listed 



in Table |3 The inverse of the mean wavenumber of the a Cet 
observations was 2.187 /im. 

As a first characterisation of the angular diameter of a Cet, 
we calculated best-fitting angular diameters ©ud and ©fdd for 
models of a uniform disc (UD: 1=1 for < fi < 1; else 0) 
and a fully darkened disc (FDD: / = //). Here, // = cos 9 is the 
cosine of the angle between the line of sight and the normal of 
the surface element of the star. The angular diameter (0ud or 
©fdd) is the only free parameter. Synthetic squared visibility 
amplitudes were calculated using the broadband VINCI sensi- 
tivity function (cf. the description in Wittkowski et al. 2004 
2006). We obtained ©ud=1 1-95 mas + 0.06 mas with a reduced 
^fy(UD) value of 3.35, and ©FDD=13.32mas + 0.12mas with 
^(FDD)= 10.20. 

Furthermore, we used an empirical parametrisation of 
the limb-darkened stellar intensity profile as I(fi) = /i a 
(cf. Hestroffer et al. I1997> . and simultaneously deter- 
mined the best-fitting values for the limb-darkened an- 
gular diameter ©ld and the limb-darkening parameter a. 
We obtain © LD = 12.27 mas + 0.05 mas and a=0.24+0.03 with 
^(LD)=1.01. This parametrisation of a limb-darkened inten- 
sity profile closely resembles the CLVs (centre-to-limb vari- 
ations) that are obtained by stellar model atmospheres based 
on plane-parallel geometry, where the atmosphere is optically 
thick for any viewing angle, and the intensity steeply drops 
to zero at the stellar limb. In this case, the Rosseland angu- 
lar diameter and our obtained limb-darkened angular diameter, 
which corresponds to the 0% intensity level, can be considered 
equivalent (cf. the discussion in Wittkowski et al. 2004 1. 

These results already show that our interferometric data 
cannot be well described by uniform or fully darkened disc 
models, but are well consistent with a limb-darkened intensity 
profile that has a strength of the limb-darkening (a=0.24) that is 
closer to that of a uniform disc {a=0) than a fully-darkened disc 
(ar=l) model. Figure^shows our measured visibility points to- 
gether with the described UD, FDD, and LD models, as well 
as with the predictions by PHOENIX and ATLAS 9 model atmo- 
spheres (derived and described below). 

3. UVES observations and data reduction 

We obtained high-resolution spectroscopy of a Cet with the 
echelle spectrograph UVES mounted on the UT2 (KUEYEN) 
telescope of the ESO VLT in service mode on 1 1 August 2003. 
We used the dichroic settings 346+580 nm and 437+860 nm 
with slit widths of 0.4" (blue arm) and 0.3" (red arm). This 
corresponds to a spectral resolution of 80 000 in the blue and 
1 10000 in the red, respectively. The details of our UVES data 
with numbers of individual integrations for each central wave- 
length of the UVES grating and exposure times are listed in 
Table|3] 

For all the observations, the bias and inter-order back- 
ground were subtracted. The spectral orders were extracted, 
flat-fielded, and wavelength calibrated with recipes avail- 
able from the ESO UVES pipeline. The "average extraction 
method" of the pipeline has been used, which uses a uniform 
average of the pixel values across the slit. Statistical errors have 
been calculated from the variance obtained with this extraction 



2003-10 






-03 4:04:48 


22.79 


67.9 


-03 4:12:01 


23.57 


68.6 


-03 4:44:19 


26.86 


70.7 


-03 4:51:27 


27.53 


71.1 


-03 5:25:02 


30.37 


72.4 


-03 5:32:18 


30.91 


72.5 


-20 6:23:10 


35.47 


72.8 


-20 6:30:14 


35.48 


72.7 


-21 5:58:25 


35.27 


73.1 


-21 6:05:40 


35.37 


73.0 


-22 3:29:06 


26.81 


70.7 


-22 3:36:30 


27.51 


71.1 


-22 4:15:05 


30.73 


72.5 


-22 4:22:10 


31.24 


72.7 


-24 5:34:53 


35.04 


73.2 


-24 5:48:19 


35.30 


73.1 


-29 5:47:00 


35.47 


72.9 


-29 5:58:10 


35.47 


72.7 


2003-11 






-27 3:14:20 


34.87 


73.3 


-27 3:21:37 


35.05 


73.2 


2003-12 






-04 1:37:03 


31.51 


72.7 


-04 1:44:08 


31.98 


72.9 


-04 2:17:53 


33.82 


73.2 


-04 2:24:54 


34.12 


73.3 


-04 3:03:23 


35.24 


73.1 


-04 3:10:23 


35.34 


73.1 


-04 4:19:06 


34.78 


71.5 


-04 4:26:25 


34.56 


71.3 


-05 0:05:53 


23.74 


68.7 


-05 0:13:52 


24.59 


69.3 


-05 1:35:55 


31.70 


72.8 


-05 1:43:08 


32.16 


72.9 


-05 2:28:03 


34.39 


73.3 


-05 2:35:39 


34.65 


73.3 


-05 3:23:51 


35.48 


72.8 


-05 3:49:52 


35.32 


72.3 


-06 3:48:38 


35.29 


72.2 


-06 3:55:49 


35.16 


72.0 


2004-01 






-12 1:00:54 


141.86 


72.7 


-12 1:13:28 


141.57 


72.4 


-13 1:28:57 


140.38 


71.9 


-13 1:40:23 


139.25 


71.6 



8.550e-01 


1.090e-02 


488 


8.322e-01 


1.073e-02 


473 


7.782e-01 


1.002e-02 


486 


7.629e-01 


9.888e-03 


467 


7.150e-01 


9.299e-03 


490 


7.052e-01 


9.208e-03 


482 


6.202e-01 


2.874e-02 


445 


6.321e-01 


2.951e-02 


414 


6.128e-01 


1.737e-02 


455 


6.062e-01 


1.718e-02 


471 


7.951e-01 


1.044e-02 


442 


7.712e-01 


9.284e-03 


425 


7.169e-01 


1.117e-02 


408 


7.073e-01 


1.302e-02 


354 


6.278e-01 


2.122e-02 


479 


6.222e-01 


2.146e-02 


351 


6.420e-01 


1.237e-02 


348 


6.573e-01 


8.384e-03 


467 


6.435e-01 


6.270e-03 


492 


6.393e-01 


6.255e-03 


488 


7.067e-01 


7.754e-03 


488 


7.033e-01 


7.740e-03 


479 


6.689e-01 


7.695e-03 


482 


6.560e-01 


7.413e-03 


476 


6.384e-01 


7.146e-03 


483 


6.367e-01 


7.200e-03 


475 


6.393e-01 


7.430e-03 


437 


6.437e-01 


9.154e-03 


388 


8.474e-01 


1.024e-02 


338 


8.037e-01 


2.597e-02 


93 


7.091e-01 


7.531e-03 


471 


6.969e-01 


7.743e-03 


467 


6.520e-01 


7.566e-03 


447 


6.548e-01 


9.280e-03 


358 


6.304e-01 


7.116e-03 


469 


6.388e-01 


7.630e-03 


424 


6.360e-01 


8.405e-03 


482 


6.376e-01 


8.533e-03 


483 


1.281e-02 


6.655e-04 


275 


1.300e-02 


6.133e-04 


299 


1.324e-02 


6.934e-04 


275 


1.365e-02 


6.689e-04 


270 
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Fig. 1. Squared visibility amplitudes and error bars of a Cet ob- 
tained with VLTI/VINCI, together with best-fitting models of 
a uniform disc (upper dashed light blue line), a fully darkened 
disc (lower dashed light blue line), a parametrisation / = fi a 
with a=0.24 (dashed blue line) and of PHOENIX and ATLAS 9 
model atmosphere predictions (solid red and green lines). The 
PHOENIX model shown has parameters r e jf=3800K, log g- 1.0, 
M = 2.3M ; the ATLAS 9 model r eff =3800K, log g= 1.0 (see 
text for more details). The upper panel shows the full range of 
the visibility function, while the middle panel is an enlarge- 
ment of the obtained squared visibility values in the first lobe, 
and the bottom panel shows an enlargement of the low squared 
visibility function in the 2nd lobe. Our measurements are sig- 



Table 3. Overview of our UVES observations of a Cet obtained 
on 1 1 August 2003 between 9:30h and 9:40h UT. Shown are the 
central wavelength /l ce ntiai of each of the UVES gratings and the 
recorded wavelength ranges A m \ n - A max . The spectra of the red 
arm of UVES were recorded using two different detectors, and 
are thus split into the two wavelength ranges. Also given are, 
for each grating, the spectral resolution R = A/AA, the number 
of individual exposures and their exposure times, as well as the 
maximum signal-to-noise ratio S /N max reached. 



■^central 

(nm) 


^min -^max 

(nm) 


R 


#of 
exp. 


Exp. time of 
each exp. (s) 


S/N nr , x 


346 


305-387 


80 000 


2 


120.0 


500 


437 


374-499 


80 000 


2 


3.7 


280 


580 


477-577/ 
583-683 


110 000 


7 


1.5 


800 


860 


665-853/ 
867-1040 


110 000 


2 


0.5 


280 



method. The extracted spectra have been flux-calibrated using 
averaged response curves as provided by the Quality Control 
Group of ESO Garching. This calibration provides a relative 
flux calibration with an accuracy of ~ 10-20% within each spec- 
trum. Absolute flux values for our UVES spectrum of a Cet 
have not been obtained. Finally, multiple exposures of each 
setting have been averaged. The resulting signal-to-noise ra- 
tios S IN vary across the spectra, and the maximum S /N values 
reached for each grating are listed in Table[5] The correction to 
heliocentric velocity has been determined by the pipeline anal- 
ysis to Vheik>=-28. 7 km/sec, so that the relation between arriving 
wavelength from the star Aq and observed wavelength A^ is 

/lobs - A (l + Vhelio/c). 

4. Atmosphere models for <i Cet 

We use new, fully line-blanketed spherical, hydrostatic atmo- 
sphere models with solar photospheric abundances (Grevesse 
& Noels [1993) obtained with version 13 of the PHOENIX code 
(Hauschildt & Baron fT999l Hauschildt et al. fT999l . Feast et al. 
( 1990) discussed that solar metallicity is appropriate for giants 
in the solar neighbourhood and found that the (rms) scatter of 
the metallicity of local giants is less than 0.08. This is con- 
sistent with the [Fe/H]= 0.00 + 0.30 derived by Decin et al. 
( 2003 ) for a Cet. The microturbulence for all our new models 
is 2kms~ 1 . The three most important input parameters for our 
spherical model are the effective temperature, the surface grav- 
ity, and the stellar mass. For each model used, we tabulate the 
flux integrated over the stellar disc from 300 nm to 1050 nm in 
steps of 0.001 nm (for comparison to our high spectral reso- 
lution UVES data). Furthermore, we tabulate monochromatic 
intensity profiles at 64 viewing angles for wavelengths from 
1.8/vm to 2.5 fim in steps of 0.5 nm (for comparison to our 
VINCI interferometric data). 

For comparison, we use intensity profiles predicted by stan- 
dard plane-parallel hydrostatic ATLAS 9 model atmospheres 
from the Kurucz CD-ROMs (Kurucz 1993 ), as in Wittkowski et 
al. J2001ll2004ll2006l . as well. The Kurucz CD-ROMs include 
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tabulated monochromatic intensity profiles for 17 angles in 
1221 frequency intervals ranging from 9.09 nm to 160 000 nm. 
In the range of the VINCI near-infrared Zf-band filter, the fre- 
quencies are sampled in steps corresponding to 10 nm. These 
data values are available as grids of effective temperature and 
surface gravity (the mass is not an input parameter for a plane- 
parallel model), and for different chemical abundances and mi- 
croturbulent velocities. For comparison, we use the grid with 
solar chemical abundances and a microturbulence of 2 km s _1 . 

For a further description of the use of the PHOENIX and 
ATLAS models for comparison to interferometric data obtained 
with VINCI, for the calculation of synthetic visibility values, 
and for effects of plane-parallel and spherical geometries on 
the obtained angular diameter, we refer to Wittkowski et al. 
120041 12006I . 

4.1. The bolometric flux of a Cet 

We have constructed a composite spectral energy distribution 
for a Cet by combining absolute spectrophotometry in the op- 
tical (Glushneva et al. 1998a 1998b) and infrared (Cohen et 
al. I1996> to estimate the bolometric energy flux, /boi, at the 
Earth (see Fig.|2}. Note that an absolute flux calibration of our 
UVES echelle spectrum has not been obtained (see Sect. [3}, 
so that we rely on spectrophotometry from other sources to 
derive the bolometric flux of a- Cet. We directly integrate the 
tabulated irradiance from 322.5 nm to 1037.5 nm (5 nm reso- 
lution) and 1.26 /mi to 35 /mi (0.02 /mi to 0.1 8 /mi resolution), 
using a five-point Newton-Cotes integration formula. Between 
1.04 /mi and 1.26 /mi, the irradiance at 20 nm intervals is esti- 
mated by linearly interpolating, in log A-log space, between 
the end points of the two data sets. We assume all the irradi- 
ance data are uncertain at the level of 5%. This procedure re- 
sults in / bo i=(1.01 + 0.05) 10~ 12 W cnT 2 . At a distance of 67 pc, 
the interstellar foreground extinction in the direction of a Cet 
is likely to be quite low. HD 18883 (B7 V), at a distance of 
134 pc, has a measured colour excess, E(B - V)=0.03, from 
space-based UV photometry (Savage et al. 1985 1 and is 0.3° 
from a Cet. Therefore, adopting E(B - V)=0.015 and an aver- 
age Ry-3.1 Galactic reddening curve (Cardelli et al. I1989i as 
a best estimate for the extinction towards a Cet, the bolometric 
flux becomes / bol =(1.03 + 0.07) 10~ 12 WcirT 2 , including the 
uncertainty in the extinction. This value is well consistent with 
the estimate by Dyck et al. ©98 1 of ( 1 .05 +0. 1 6) 1 0" 12 W cnT 2 
based on broadband photometry alone. 

4.2. Comparison of PHOENIX model predictions to 
available spectrophotometry 

We have compared synthetic spectral energy distributions 
(SEDs) from PHOENIX models with the absolute spectropho- 
tometry used to derive the bolometric flux above. The goodness 
of the SED fit is largely insensitive to specific values of the 
surface gravity and mass in the model grid. The least-squares 
fits for the grid and application of the F-test provide the fol- 
lowing mean values and 1 <x uncertainties: r e ff=3800+ 100 K 
for 0r oss =12.O ± 0.7 mas. These uncertainties are correlated, 



so that low r e ff values correspond to high Qr oss values and 
vice versa. Figure |2] shows the spectrophotometric data com- 
pared to the synthetic SED based on the model with mean 
parameters mentioned above. Differences between observed 
and synthetic SED are largest at optical (A < 0.7 nm) wave- 
lengths. With the Hipparcos parallax, the angular diameter of 
®Ross=12.0 + 0.7 mas corresponds to a Rosseland linear radius 
of # Ross =87 + 10/?o and luminosity log(L)=3.16 + 0.18L o . 
Relying on the evolutionary tracks from Girardi et al. (J2000 1, 
the position of the star in the theoretical HR diagram results in 
a mass M/M Q =2.3 + 0.5 (see Fig. 1 in Wittkowski et al. l2004> . 
From the mass and the Rosseland radius, we derive a surface 
gravity of log g-0.9 + 0.2. 

4.3. Comparison of model predictions to our VINCI 
data 

We derive best-fitting limb-darkened (0% intensity) diame- 
ters ©ld to our VINCI data for different PHOENIX model at- 
mospheres with parameters T e g, logg, and M that lie within 
ranges consistent with the estimates above. The procedure used 
is the same as described in Wittkowski et al. (2004 2006 1. 
Table0]lists the considered model parameters and the resulting 
best-fitting ©ld, ©Ross, and reduced^ 2 values for each model. 
The ratio between the outermost model layer corresponding 
to ©ld and the Tr oss =1 layer corresponding to ®r oss is read 
from the PHOENIX model atmospheric structure as described in 
Wittkowski et al. (2004 1. The precision of the best-fitting di- 
ameter values for a given model atmosphere listed in Table@]is 
relatively high (0.4%). We have shown by a study of night-to- 
night variations that such high-precision diameter values based 
on VINCI data are generally reliable (Wittkowski et al. 2006 1. 
Also, the different considered PHOENIX models with variations 
of T e g , log g, and M result in consistent Rosseland angular di- 
ameters, so that no uncertainty due to the choice of model pa- 
rameters needs to be added. 

For comparison, we derive the best-fitting limb-darkened 
diameter based on the corresponding plane-parallel ATLAS 9 
model as well. For plane-parallel geometry, ©ld and @r oss can 
be considered equivalent (Wittkowski et al. 2004 1. However, 
the definition of the Rosseland diameter is more precise for 
models based on spherical geometry (cf. Wittkowski et al. 
2006 1. Here, the resulting Rosseland angular diameters based 
on the PHOENIX and ATLAS 9 models agree very well. 

The corresponding synthetic /T-band visibility curve and 
thus the resulting xl values are not sensitive to variations of 
r e ff, logg, and M within the considered ranges. The good xl 
values of 0.99 shows that the measured and model-predicted 
strength of the limb-darkening effect in the /f-band agree well. 
Figure [2 shows the measured squared visibility amplitudes to- 
gether with the PHOENIX and ATLAS model prediction. Also 
shown are the simple models of a uniform disc, fully darkened 
disc, and the parametrisation / = /j" as discussed in Sect. [2 

The high-precision Rosseland angular diameter obtained 
from the comparison of the PHOENIX model atmospheres 
to our VINCI data of Qr oss =12.20 + 0.04 mas corresponds 
with the Hipparcos parallax to a linear Rosseland radius of 
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Fig. 2. (Top) The observed spectral energy distribution of a Cet (shown as error bars) from Glushneva et al. ( 1998a, 1998b) and 
Cohen et al. {1996 1 compared with synthetic spectrophotometry ('x' symbols) derived from the PHOENIX model (shown in grey 
at high-resolution) with the mean parameters listed in the text. (Bottom) The percentage difference between the observed and 
synthetic spectrophotometry in each wavelength bin. 



Table 4. Best fitting angular diameters of a Cet based on a com- 
parison of our VINCI data to different model atmospheres. 



Model parameters 




©LD 




©Ross 


A v 


T eB [K] 


l°g£ 


M/M 




[mas] 




[mas] 




PHOENIX models: 














3800 


1.0 


2.3 


12.60 ± 0.04 


12.20 


±0.04 


0.99 


3700 


1.0 


2.3 


12.60 


±0.04 


12.20 


±0.04 


0.98 


3900 


1.0 


2.3 


12.60 


+ 0.04 


12.20 


±0.04 


0.99 


3800 


0.5 


2.3 


12.95 


±0.04 


12.19 


±0.04 


0.99 


3800 


1.5 


2.3 


12.43 


±0.04 


12.21 


±0.04 


0.99 


3800 


1.0 


2.0 


12.63 


±0.04 


12.20 


±0.04 


0.99 


3800 


1.0 


2.6 


12.58 


±0.04 


12.20 ± 0.04 


0.99 


ATLAS 9 model for compari 


son: 










3750 


1.0 


/ 


12.20 


±0.04 


12.20 


±0.04 


0.99 



^Ross=89 ± 5 R Q . Here, the error is by far dominated by the er- 
ror of the bolometric flux. Our VINCI value for 0r oss together 
with the bolometric flux derived above results in an effective 
temperature of r e jf =3795 + 70 K. /?r oss and r e fr result in a lumi- 
nosity logL/L =3.16 + 0.08 (L=146OL ± 300 L ). These val- 
ues are well consistent but more precise than the values es- 
timated above by a comparison of the PHOENIX models with 
available spectrophotometry. Based on the stellar evolution- 
ary tracks by Girardi et al. (2000), we can estimate a mass of 
MIM Q -23 ± 0.2 and a surface gravity of log g-0.9 + 0.1 (cgs). 



4.4. Comparison of model predictions to our LIVES 
spectra 

We consider in the following the spherical PHOENIX model 
atmosphere as described above with parameters r e ff=3800K, 
logg=1.0, and M = 2.3M because it is closest to the model 
parameters derived by comparison to our VINCI data and well 
describes the measured strength of the limb-darkening in the 
/T-band. Here, we compare the model-predicted spectrum in 
several bandpasses to our measured UVES spectrum to inves- 
tigate to what extent the same model that is consistent with 
our interferometric data is consistent with our high-resolution 
spectroscopic data as well. 

We use four dominant TiO bandheads around 5598A, 
7054A, 7088A, and 7126Afor our analysis. In addition, we 
have selected 6 bands around dominant atomic lines of Fe I 
(3683 A, 5447 A, 6945 A), Cal (4227 A), and Til (5966A, 
9675A) that cover all gratings used for our UVES spectrum and 
span a total range from 360 nm to lOOOnm. Profiles of atomic 
lines and TiO bandheads such as these are good indicators of 
stellar model atmosphere parameters, most importantly, effec- 
tive temperature and surface gravity (e.g., Valenti et al. [1998 1. 



Calibration of spectral resolution and wavelength The wave- 
length scales of the synthetic spectra were corrected for the air 
refraction index, for the Earth's motion and rotation, and for 
the radial velocity of a Cet of v ra d=-26. 08 km/sec (Famaey et 
al. 2005 ). For each bandpass separately, the synthetic PHOENIX 
spectrum was broadened by convolution with a rotational pro- 
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Fig. 3. Observed UVES spectrum (black) together with the synthetic spectrum (red) for each bandpass considered. The parameters 
of the spherical PHOENIX model used are r e ff =3800 K, log g- 1 .0, M = 2.3 M Q , and solar chemical abundance. The vertical dotted 
lines denote the central wavelength and bandpass used to compute the^J and integral W b s /W mo d e i values in Table|5] The vertical 
dashed lines denote the bandpasses blueward and redward of the central bandpass that are used to adjust the scale factor / 
between model spectrum and measured spectrum and to normalise the equivalent width. 



8 



M. Wittkowski et al.: Interferometry and spectroscopy of Menkar 



Table 5. Bandpasses used for the comparison of PHOENIX 
model predictions to our UVES spectrum. Listed are the cen- 
tral wavelength /l a ; r , the identification (ID), the obtained wave- 
length shift AA, the scale factor / between observed and model- 
predicted spectrum, the reduced xl value, and the ratio of ob- 
served and model-predicted equivalent width. The scale fac- 
tors / are derived in blue and red bandpasses at [-8.. - 4] A 
and [4. .8] A with respect to the central wavelength. The^J val- 
ues are then derived in bandpasses [-2..2]A around /l a ; r . The 
equivalent widths are the integrals of the central bandpasses af- 
ter normalisation of the spectrum to unity using the integral of 
the blue and red bandpasses. 





ID 


AA 


./norm 


X v 


W^obs / Wmodel 


A 




A 








5598.41 


TiO 


-0.031 


0.89 


13378 


0.70 


7054.19 


TiO 


-0.006 


1.09 


530 


1.01 


7087.60 


TiO 


-0.010 


1.04 


1071 


0.80 


7125.59 


TiO 


-0.029 


1.07 


674 


0.91 


3683.06 


Fel 


-0.002 


1.09 


17496 


3.4 


4226.73 


Cal 


-0.004 


0.68 


169 


0.91 


5446.92 


Fel 


-0.009 


0.87 


32299 


0.47 


5965.83 


Til 


-0.037 


0.84 


7051 


0.96 


6945.20 


Fel 


-0.007 


1.20 


718 


0.69 


9675.49 


Til 


-0.073 


1.21 


188 


0.79 



file using v sin i-3 km/sec, which is the mean value for M0 
giants (Munari et al. 2001 1. It was also broadened with a 
Gaussian profile to match the spectral resolution of the obser- 
vation given in Table[3] Due to possible residuals of aCet's ra- 
dial velocity, and residuals of the synthetic line positions, small 
offsets between the synthetic and observed wavelength scales 
can remain. We have used the IDL routine crscor from the 
IUE library (http://archive.stsci.edu/iuei to cross-correlate the 
observed and synthetic spectra for each bandpass used and to 
derive remaining wavelength shifts. The resulting wavelength 
shifts as listed in Tableware of the order of 0.01 A (correspond- 
ing to 0.6 km/sec at 5000A) and are small for the purpose of our 
comparison. 

Comparison of observed and synthetic spectra An absolute 
flux calibration of our UVES spectra was not obtained, lead- 
ing to an unknown scale factor between observed and synthetic 
spectrum. We have determined, by a standard least-squares fit, 
the best scale factor / to match observed and synthetic spec- 
tra in two bandpasses blueward and redward of each central 
bandpass. For a cool giant, we cannot find good continuum 
bands close to the central wavelengths. The bandpasses used 
for this adjustment have a standard width of 4A and standard 
central wavelengths -6A and +6A with respect to the central 
wavelength, and inevitably contain spectral features as well. 
The resulting xl value between observed and synthetic spec- 
tra was calculated for each central bandpass. The normalised 
factors /norm and the^J values are listed in Table [5] The vari- 
ations of /norm are generally consistent with the relative uncer- 
tainty of the flux calibration of our UVES spectrum of 10-20% 
for each spectrum as quoted in Sect. [3] In addition to the di- 



rect x 2 comparison between observed and synthetic spectra, 
we calculate the ratio between measured and predicted equiva- 
lent width W b s /W moc iei as a measure of the consistency of the 
integrated line strengths. The equivalent width is derived, sep- 
arately for synthetic and observed spectra, by integrating the 
spectrum over the central bandpass after normalisation to unity 
using the integral of the red and blue bandpasses. 

Figure |3 shows the final comparison of observed and syn- 
thetic spectra for each considered bandpass. Our PHOENIX 
model describes the measured spectrum around all selected fea- 
tures well qualitatively. We have inspected a number of addi- 
tional features across the total wavelength range of our UVES 
spectrum and find comparable agreement. However, on the de- 
tailed level of the high spectral resolution obtained, measure- 
ment and model prediction exhibit differences that dominate 
the resulting xl value in Table[5] These differences include dif- 
ferent strengths of individual lines, line positions, and features 
that appear in the observed spectrum, but not in the synthetic 
spectrum, and vice versa. It is known that synthetic and ob- 
served stellar spectra at this detailed level show discrepancies 
such as these and that relatively large xl values are thus to be 
expected for a direct comparison (see, e.g., Kucinskas et al. 
2005). The integrated quantity W b s /W mo d e i has a mean value 
of 0.80 and standard deviation 0.16 (excluding the 3683.06 A 
Fe I bandpass where the blending of different spectral features 
is strongest), i.e., the observed spectral features have, over all, 
a lower strength than predicted by the model at the level of 

1.25 cr. 

A variation of model parameters log g, and M within the 
uncertainties that remain from the comparison to our VINCI 
data do not lead to significant changes of Fig. [3] and on the 
xl and Wobs/Wmodei values in Table [3] These parameters can 
thus not be further constrained by the comparison to our 
UVES spectrum. A variation of T e g (model values 3750 K, 
3800 K, 3850 K, 3900 K; other model parameters unchanged) 
allows us to estimate the effective temperature of aCet to 
r e ff ^3820 + 50 K both based on the or the Wobs/Wmodei val- 
ues as a function of T e g. This estimate is based on the com- 
parison of observed and synthetic spectral lines/bands and is 
consistent with the estimate based on the best-fitting diameter 
derived from the interferometric data and the bolometric flux 
derived from available spectrophotometry as described above. 

There are three major limitations to a stronger constraint 
of the model parameters. Firstly, an absolute flux calibration of 
the UVES spectrum was not obtained, and the relative flux cal- 
ibration within each spectrum reaches a precision of not better 
than 10-20%. Secondly, for cool giants, the selected spectral 
features are not isolated and the bandpasses used are inevitably 
contaminated by several other lines; nearby true continuum 
bandpasses are not available. Finally, the total differences be- 
tween observed and synthetic spectrum for our bandpasses, as 
characterised by the^J values in Tabled are dominated by de- 
tailed effects other than the main model parameters T e g , log g, 
and M. 
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Fig. 4. Correction factors from the 0% LD diameter Old to 
the UD diameter based on our derived PHOENIX model for 
a Cet with parameters r eff =3800K, logg=l, and M/M =2.3. 
The dashed line denotes the Rosseland diameter 0r oss with 
respect to the 0% intensity diameter ©ld, which is 0.968 for 
this model. The thin grey line shows the full resolution of our 
PHOENIX model file (5 A), and the solid black line is reduced to 
a resolution of 20 nm. 

Table 6. UD diameters predicted for several previously used 
bandpasses based on our derived PHOENIX model for a Cet 
with parameters r e ff=3800K, logg=l, M/M Q -2.3, solar chem- 
ical abundance, and 0LD=12.6Omas. We estimate the accuracy 
of the model predictions to be < 5%, see text for details on the 
calculation. 



Ac 


6A 


W UD 




W UD 


Reference 


nm 


nm 


mas 




mas 




2187 


400 


11.90 


11.95 ±0.06 


Direct fit in this work 


450 


20 


11.24 


11.33 ± 


0.29 


Mozurkewich et al. I1991i 


451 


20 


11.24 


11.33 ± 


0.41 


Mozurkewich et al. 1 2003 1 


550 


20 


11.09 


11.47 ± 


0.25 


Mozurkewich et al. 1 2003 1 


712 


12 


11.31 


11.95 ± 


0.23 


Quirrenbach et al. 1 19931 


754 


5 


11.52 


11.66 ± 


0.22 


Quirrenbach et al. 119931 


800 


20 


11.60 


12.25 ± 


0.16 


Mozurkewich et al. 1 1991 1 


800 


20 


11.60 


12.27 ± 


0.24 


Mozurkewich et al. 120031 


2200 


400 


11.90 


11.6 


t 0.4 


Dyck et al. fPMl 



Table 7. Fundamental parameters of the M 1 .5 giant a Cet 
based on the analysis in this paper. The mass and surface grav- 
ity rely on the evolutionary tracks by Girardi et al. (2000 1. 



Parameter 



Value 



Rosseland angular diameter 
Rosseland linear radius 
Bolometric flux 
Effective temperature 
Luminosity 
Mass 

Surface gravity 



© Ross = 12.20 ±0.04 mas 

fl Ross = 89 ± 5R Q 

/ b0 , = (1.03 ± 0.07) x 10- 12 W/m 2 

T cff = 3795 ±70 K 

logL/L = 3.16 ±0.08 

M = 2.3 ± 0.2M o 

logg = 0.9±0.1 (cgs) 



ibility data were obtained at other parts of the visibility curve. 
We average the correction factors over rectangular bandpasses 
with central wavelength A c and width 6A. The exact shape of 
the bandpasses used for the previous observations and the exact 
signal processing of the different instruments is not taken into 
account. We estimate the total error resulting from these sim- 
plifications to be less than about 5%. Figure|4]shows the result- 
ing correction factors for the wavelength range from 400 nm 
to 2500 nm. Table [6] lists the previously obtained interfero- 
metric UD diameter measurements of a Cet compared to the 
prediction by our derived PHOENIX model for each specific 
bandpass used. Only the measurements at 550 nm and 800 nm 
by Mozurkewich et al. ( 1199 II 12003 ) and the measurement at 



4.5. Comparison of model prediction to previous 
diameter measurements 

Several UD diameter measurements of a- Cet based on long- 
baseline interferometry have previously been obtained (see 
Sect. Here, we use our derived PHOENIX model of a Cet 
(Teff=3800K, logg=l, M/Mq-23, solar chemical abundance, 
©LD=12.60mas) to predict the UD diameter at the previously 
used bandpasses. The correction factor from the 0% (LD) di- 
ameter to the UD diameter is calculated so that both visibil- 
ity curves, the UD curve and the curve based on the model 
atmosphere, match where |V| 2 = 0.3, as previously done by, 
e.g., Hanbury Brown et al. ( 1974 1, Quirrenbach et al. (0996 1, 
and Wittkowski et al. (J2001 1. The choice to match the visibil- 
ity curves at |V| 2 = 0.3 is arbitrary, and the accuracy of the 
resulting correction factor may decrease if most measured vis- 



7 12 nm (TiO band) by Quirrenbach et al. ( 1993 1 are not consis- 
tent within 1 <x with the nominal value of our model prediction, 
but still within the estimated ~5% uncertainty of our predic- 
tion. The discrepancy at the TiO bandpass may either be caused 
by an imperfect match of the exact instrumental bandpass and 
signal processing (as here the strength of the limb-darkening ef- 
fect changes rapidly), or by an imperfect model description of 
the spatial structure of the layers where TiO molecules reside. 
It has been reported that the use of different line list combina- 
tions of molecules lead to significantly different model struc- 
tures and spectra, in particular in the optical where TiO bands 
are important (Allard et al. 2000, Kucinskas et al. 2005). 



5. Summary and discussion 

We have shown that our derived PHOENIX model atmo- 
sphere for a Cet (Menkar) with model parameters r e ff =3800 K, 
logg=l, and M/M =2.3 is consistent with both the measured 
strength of the limb-darkening in the broad near-infrared K- 
band and with the profiles of spectral bands around selected 
atomic lines and TiO bandheads from 370 nm to 1000 nm. 
At the detailed level of our high spectral resolution (R up to 
1 10000), however, noticeable differences of observed and syn- 
thetic spectra exist. The discrepancies include differences in 
the strengths and positions of spectral lines/bands, as well as 
detailed spectral features that only appear in either the ob- 
served or synthetic spectrum. It has previously been reported 
that the existence of such detailed effects, especially in the op- 
tical where TiO bands are important, is known and may be due 
to, for instance, different line list combinations or treatments of 
convection (Allard et al. 2000 Kucinskas et al. 2005). 
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Fig. 5. Measured effective temperature of the M giants ySge 
(Wittkowski et al. 2006 1, cCet (present paper), and iftPhe 
(Wittkowski et al. 2004 1 versus spectral type compared to em- 
pirical calibrations for giants. The three measurements are con- 
sistently obtained by Rosseland diameters based on a compari- 
son of VINCI interferometric data to spherical PHOENIX model 
atmospheres and bolometric fluxes obtained by careful integra- 
tion of available spectro-photometry. 



Our results illustrate as well the power of combining in- 
terferometry and high-resolution spectroscopy to constrain and 
calibrate stellar model atmospheres. The newly offered instru- 
ment AMBER (Petrov et al. 120031 at the VLTI permits the 
recording of spectro-interferometric data with a spectral res- 
olution of up to 12000, i.e., with a resolution that is several 
orders of magnitudes above that of previous interferometric in- 
struments. This spectral resolution of the newest generation of 
interferometric instruments, however, is still clearly below that 
of modern optical spectrographs, such as UVES, or IR spectro- 
graphs, such as CRIRES, which provide a spectral resolution 
of up to R > 100000. It will thus be very valuable to combine 
AMBER interferometry with UVES or CRIRES highest reso- 
lution spectroscopy to constrain stellar models. A limitation of 
our current study is the lack of an absolute flux calibration of 
our UVES data, which should be obtained with high precision 
for such suggested future studies. 
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The fundamental stellar parameters of a Cet are most con- 
strained by our high-precision angular diameter obtained from 
the comparison of the PHOENIX atmosphere model with our 
VINCI data, and the bolometric flux based on available spectro- 
photometry. This resulting set of fundamental stellar parame- 
ters of a Cet is summarised in Table These values are con- 
sistent with previous measurements mentioned in Sects. H an d 
14.51 but generally have a higher precision. 

Together with the results on the M4 giant i/rPhe 
(Wittkowski et al. 2004] Paper II of this series), the M0 gi- 
ant y Sge (Wittkowski et al. 2006, Paper III of this series), and 
the M 1 .5 giant a Cet (present paper), we have consistently ob- 
tained sets of high-precision fundamental parameters for these 
three M giants based on a comparison of VINCI data to spher- 
ical PHOENIX model atmospheres. Figure[5]shows the obtained 
effective temperatures versus spectral types compared to em- 
pirical calibrations by Ridgway et al. ( 1980), Di Benedetto & 
Rabbia ( 1987), and Dyck et al. ( 1996). The measurements are 
consistent with the empirical calibrations shown at about the 
1 <t level, but may suggest a flatter temperature slope for early 
M giants and a stronger decrease of effective temperature for 
cooler (toward M 4-5) giants. 

The studies on these three cool giants also consistently 
show that the model-predicted strength of the limb-darkening 
is not significantly affected by the geometry of the model at- 
mosphere for the broad near-infrared Zf-band, and also for nar- 
rower visual bandpasses in the case of the NPOI observations 
of ySagittae of Paper III, with the precision of current mea- 
surements. Effects from plane-parallel versus spherical geom- 
etry thus appear to be more subtle than generally seems to be 
expected. However, the most precise definition of a meaning- 
ful diameter, such as the diameter where the Rosseland-mean 
optical depth reaches unity, can best be obtained based on a 
spherical model, as discussed in the articles of this series. 
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